% 创建时间：2025-06-20 23:10:39
% 作者：兪馳
clc;
clear;
A=[1 2 3 4;
    2 2 1 1;
    2 4 6 8;
    4 4 2 2]; 
B=[1:4]'; 
C=[A B]; 
rank(A), 
rank(C)
x=pinv(A)*B, norm(A*x-B)